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Abstract 

We study low-energy virtual Compton scattering off the proton within 
the framework of a nonrelativistic constituent quark model. The Compton 
tensor is divided into two separately gauge-invariant contributions. The first 
consists of the groundstate propagation in the direct and crossed channels 
together with an appropriately chosen term to satisfy gauge invariance. The 
residual part contains the relevant structure information characterized by 
the so-called generalized polarizabilities. We discuss two different schemes 
to obtain the generalized polarizabilities from the residual term. Explicit 
predictions for the generalized polarizabilities are presented for the Isgur- 
Karl model. Our results are compared with previous predictions in that 
model as well as other approaches. 
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Typeset using REVT^ 



I. INTRODUCTION 



The electromagnetic polarizabilities a (electric) and (3 (magnetic) of real Compton 
scattering (RCS) [I]] describe, to leading order in the photon frequency, the model- 
dependent response of a spin-0 or spin-1/2 system beyond the low-energy theorem (LET) 
l^-Q. Within a classical framework these quantities are accessible to an intuitive inter- 
pretation as a "measure of the stiffness or rigidity of a system p." There have been 
considerable experimental efforts to determine the proton polarizabilities from Compton 
scattering off the proton Until recently, the most precise values for the proton 

polarizabilities were derived in the work of MacGibbon et al. and analyzed by means 
of dispersion relations at fixed t [^]. The results were otp = (12.1 ± 0.8 ± 0.5) x 10~^ fm^ 
and Pp = (2.1 =1=0.8 =F 0.5) x 10""^ fm^. The analysis of Ref. made use of the Baldin sum 
rule |TT[| which relates the sum of the two polarizabilities to the total photoabsorption cross 
section. The sum rule constraint was ap + f3p = (14.2 ± 0.5) x 10~^ fm^ |jl2|, while more 
recent analyses found (13.69±0.14) x 10"^ fm^ [||] and [14.0± (0.3-0.5)] x 10"^ fm^ [|14 



New real-Compton-scattering data below pion-production threshold have been measured 
at the Mainz Microtron (MAMI). The new global fit including these data results in the 
values ap = (12.24±0.24±0.54±0.37) x 10"^ fm^ and pp = (1.57±0.25±0.52±0.37) x 10"^ 
fm^ |15|]. As there is no free neutron target the experimental information on the neutron 
polarizabilities is much less certain. Results for the electric polarizability have been ob- 



tained from low-energy neutron scattering off the Coulomb field of a heavy nucleus [16-19 



Alternatively, the quasifree Compton scattering reaction (i(7, '-f'n)p [pil|-p2| as well as elas- 
tic deuteron Compton scattering |23] have been investigated to extract information on the 
neutron polarizabilities. 

Clearly, the concept of polarizabilities is open to generalizations in, at least, two 
directions. One possibility consists of investigating higher-order terms in the expansion of 
the RCS amplitude. For example, for a spin-1/2 system one finds four additional constants 
(spin polarizabilities) at third order |24-p7| and four terms for the spin-averaged amplitude 
at fourth order |2^J2^. Another option is to allow at least one photon to be virtual. This 
second generalization, which was already discussed in the 1970's in the context of nuclei 
2^, has recently been applied to the nucleon pil| , pT[ . 



As in all studies with electromagnetic probes, the inclusion of virtual photons substan- 
tially increases the possibilities to investigate the structure of the target. The use of a 
virtual photon allows one to access longitudinal degrees of freedom and to vary the three- 
momentum and energy transfer to the target independently. The potential of investigating 
electron-proton bremsstrahlung as a source of the virtual Compton effect amplitude was 
already noticed by Berg and Lindner [^] in the 1960 's. In its most general form for one 
real and one virtual photon, this amplitude can be described by twelve form factors of 
three invariants [^. A generalized low-energy theorem (GLET) analogous to the LET of 
RCS was derived in Ref. where it was shown that, up to and including second 
order in the momenta, all twelve amplitudes for the proton can be expressed in terms of 
the proton mass M, the anomalous magnetic moment k, the electromagnetic Sachs form 
factors Ge and Gm, the mean square electric radius r"^, and the RCS electromagnetic 
polarizabilities a and (3. 

In Ref. [pOl, the model-dependent response beyond the LET was analyzed by means 
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of a multipole expansion. Only terms contributing to first order in the frequency of tlie 
outgoing real photon were kept, and the result was expressed in terms of ten generalized 
polarizabilities (GPs) which are functions of the three-momentum of the virtual photon 
in the initial state. Further progress has been made with respect to implementing the 
constraints due to the discrete symmetries in combination with particle crossing. In Refs. 
33-RB| it was shown that only six of the originally ten GPs are independent if charge- 
conjugation symmetry is combined with particle crossing. 

Following the very first calculation in the framework of a nonrelativistic quark model 



30|j37[| , there have been numerous predictions for the GPs within various approaches, 



including phenomenological Lagrangians |]38| , |39| , the linear sigma model |^0|5 chiral per- 
turbation theory [^, the Skyrme model a relativistic constituent quark model |^ 



and the so-called small scale expansion |Q. On the experimental side, first evidence for 



virtual Compton scattering events was reported in Ref. |^5[. The first results for two 
structure functions involving linear combinations of GPs at = 0.33 GeV^ have been 
obtained from a dedicated VCS experiment at MAMI Further experiments probing 
the GPs at different values of momentum transfer are underway at Jefferson Lab |^ and 
MIT-Bates ISl. 



Our work is organized as follows. We start out in Sec. II with a general discussion 
of the (virtual) Compton scattering tensor for a nonrelativistic system of particles. 
We propose a separation into two individually gauge-invariant pieces. The first part 
consists of the groundstate propagating in the direct and crossed channels, supplemented 
by an appropriately chosen term to satisfy gauge invariance. The residual part contains 
the structure information contributing to the GPs. Section III deals with the multipole 
expansion of the Compton tensor and the definition of the GPs according to Ref. |]30[ . 
Two schemes are presented to identify those pieces of the residual term that actually 
yield contributions to the GPs. In Sec. IV we reconsider the calculation of the GPs in 
the framework of the Isgur-Karl model, and Sec. V contains a short summary and some 
conclusions. The more technical details can be found in the appendices. 



II. HADRONIC COMPTON SCATTERING TENSOR 



In this section we discuss the formalism to describe both real and virtual Compton 
scattering 'j*{uj, q) + N{Ei,pi) 7*(i^', q') + N{Ef,pf) off a composite system, denoted 
by A^, within the framework of nonrelativistic quantum mechanics. We only consider the 
main steps to derive the hadronic tensor of Compton scattering and refer the interested 
reader to Refs. p9|-p^ and the more recent work of Refs. p^j55[] for a detailed discussion. 

The starting point is the nonrelativistic Hamiltonian for a composite system of A^ 
particles with masses rria, 



Ho 



N 

E 



2m„ 



a<l3 



where fa and —i Va refer, respectively, to the position and conjugate momentum of the 
particle a in the coordinate representation. We use natural units, i.e., h = c = 1, e > 0, 
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and a = e^/47r ^ 1/137. For simplicity, we only consider a local potential i.e., 
the interaction between the constituents does not contain momentum-dependent and/or 
exchange forces, and thus avoid the problem of exchange currents. The eigenstates of Eq. 
(0), denoted by iXpx), are normalized according to 

{X'px'lXpx) = 6''{px'-px)Sx'x, (2) 

where X and px refer to internal quantum numbers and the total momentum of the 
system, respectively. 

The interaction with an external electromagnetic field is introduced via minimal sub- 
stitution. In addition, we include a coupling to an intrinsic magnetic dipole moment 
fla = ea^af^^a of the constitucuts, where aQ,/2 is the Pauli spin operator. The resulting 
electromagnetic interaction Hamiltonian in the Schrodinger representation readsQ 

Hi = Hi^, + Hi^2, (3) 
Hj^,= I SxJ^'{x)A^{x), (4) 



where A'^{x) is the second-quantized photon field and 



N 



m = E 



(6) 

N 

P(^) = eaS^{x - fa), (7) 



0=1 



N 2 

S^o = B"^' = 0, B'^ix, x') = J2 - ra)6%x' - f„ 



a=l 



The total Hamiltonian H = Hq + Hj-ad + Hj, where Hrad refers to the Hamiltonian of the 
radiation field, is time independent. 

The hadronic Compton tensor is obtained by calculating the contributions of Hj 2 and 
Hi I in first-order and second-order perturbation theory, respect ively,^ 

Mj^^(g', q,p) = S^::{q', q) + Tf^iq', q,p), (9) 

where q and q' refer to the four-momenta of the initial and final photons. In Eq. (P) 
we have also kept the dependence on a third independent three-momentum, namely, the 



We adopt a covariant notation, even though a calculation within the framework of nonrela- 
tivistic quantum mechanics does not provide a covariant result. 

^For notational convenience we keep the four-momenta q and q' as arguments, although their 
time components are related by energy conservation. 
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average of the initial and final target momenta p = {pi + pf)/2. Usually this dependence 
is reexpressed in terms of the photon three-momenta when specifying the reference frame 
through p = {{q — g')/2,0, — (g + q')/2} for the laboratory, Breit, and center-of-mass 
frame, respectively. This issue is of importance in the context of photon crossing as will 
be shown later. 

Our normalization convention for the Compton tensor follows from the S-matrix ele- 
ment of RCS, 7(g, A) + N{pi, a) ^ 7(5', A') + N{pf, a'): 

Sj, = 5\pf - p.)5\q- q')S^>Jyx " —^6\pf + q' - p. - q)e';{\')e,{X)M^^ {q\ q, p). 

(10) 

Using the interaction representation with respect to Hq + Hrad^ the first-order con- 
tribution is a sum of the "seagull" terms, 

^;r(g',9) = {2Txf{fpj\B^^{q\q)W>), (11) 

where 

E'^^(g', q) = \j dh [e-^^'-'B^^'iz, 0) + e"^-'B''^{z, 0)] . (12) 
The latter equation follows from translational invariance, 

B^"'{x,x') = e-'^-""E'^"(x - x',0)e+'^-"", 

where P refers to the total momentum operator of the composite system. 
The second-order result from the direct and crossed channels reads 

T,7(g',g,p) = {2nfY. f d'px{fPf\Jm\Xpx) {Xpx\rm^V^) 

X -I ^fVPf) + - ^x{Px) 

I d'px{m\j'mxpx) ^^ff.^'^'', ~ ^; ^^V. {xpx\nmv.). 

X ^f\Pf) - ^ - ^xypx) 

(13) 

with the energy of the intermediate state \Xpx) given by Exipx) = Px/2M + /S.Ex, M 
being the target mass and Ai?x the excitation energy of the corresponding state. For 



further reference we note that Eqs. (11) and (|T^) are symmetrical under photon crossing, 
q^ —q'^ and u /i.R 



^Alternatively, we could have treated the electromagnetic field as a time-dependent external 
field and used the interaction representation with respect to Hq. Both approaches yield the 
same result. 

^ Sometimes the symmetry property B^^ {x,x') = B^^[x' ,x) is used to express Eq. (|l2|) as 
/ d^ze~^^ '^B^'^{z, 0) which, however, is no longer manifestly symmetrical under photon crossing. 
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By use of the standard procedure to separate center-of-mass and internal motions 
|]5T|J53[| , the matrix element of the current operator may be written as 

ifpjlUOW,) ^ ^(/l J^(g,p)|^), (14) 

where 

^o(g,p) = (15) 

with q = Pf — Pi- In Eq. ( |l5|) , the intrinsic current operator i^^{q) and the charge density 
operator p(g) are 

r{q) = E {{e''-''-,p'^] + ^Sc. X qe^'^-''A , (16) 
p(g) = Ee.e^^^-^^^ (17) 

a, 

where {a, h} = ah + ha is the standard anticommutator, and r and p'^ are the intrinsic 
coordinates and momenta of the particles relative to the center of mass, 



f'^ = fo,-R, ^=-^^mar;, (18) 

P'.=Pa-'^P, P = Y.P.. (19) 



We note that the intrinsic coordinates and momenta do not satisfy the canonical commu- 
tation relations, 

[r'.M=^8^,{5.p-'^). (20) 

From now on we assume that both initial and final states, denoted by |0) and (0'|, 
respectively, correspond to the ground state of the system. However, we explicitly allow 
for a change in the spin projection. Inserting Eq. (P^ ) into Eqs. (|TT|)-(|T^ and integrating 
over the center-of-mass momentum, the Compton tensor of Eq. (^ can be cast into the 
formQ 



X 



(21) 



e? 



^In the following we omit the subscript fi. 
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Before calculating the GPs, we split the Compton tensor into two parts that are 
separately gauge invariant and symmetrical under photon crossing, 

M^"" = M'^" + M^. (23) 

The modified pole term Mp" is defined as 

where the pole term Tp' corresponds to the contribution of the intermediate ground state 
in the direct and crossed channels and G^'^ reads 

G^\q',q) = G'''{q\q) = 0, G^^{q',q) = <5,,.-L(0'|p(-g')|0')(0|p(g )|0)5m,m.. (25) 

A derivation of G^" is given in Appendix A. In particular, Mp' generates the correct 
Thomson amplitude in the limit q, q' 0. 
The residual term Mp^ is then given by 

where Tj^'^ is the contribution of the excited states in the direct and crossed channels, 
while the modified seagull term is S^'^ = S'^" — G^" . 



III. RESIDUAL COMPTON TENSOR AND GENERALIZED 

POLARIZABILITIES 

In this section we review the multipole expansion of the residual tensor and define the 
GPs according to Ref. In this context, we generalize the discussion of Ref. by 
allowing for a general spin / of the initial and final states, respectively. We then discuss 
two different schemes of evaluating the residual term Mj^ , the first one based on the 
presentation of Ref. |^ and the second one on a systematic 1/M expansion, which is 



capable of incorporating the constraints of photon-crossing symmetry to leading order in 



A. Multipole expansion and generalized polarizabilities 



The starting point of Ref. |3y] for the definition of the GPs is the multipole decompo- 
sition of the residual term, which is to be evaluated in the center-of-mass frame defined 
by Pi + q = Pf + q' = 0,13 



M^''{Mf,uj',q';M„uJ,q) = 

47r E g,','V^{pL'M',q')Hi^''''^''''''''\Mf,u;', \n M.,u;, \q\)V''*ipLM,q)g^ 



ppi 



(27) 



p,L,M, 
p' ,L',M' 



^We repeat that uj and uj' are related by energy conservation. 
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where {V^{pLM, q)} constitutes the four-dimensional basis of the multipole expansion of 



four- vector fields as introduced in Appendix C of Ref. |3^. In particular, p(p') denotes 
the type of the initial (final) multipolarity (p = scalar, p = 1 magnetic, p = 2 electric, 
and p = 3 longitudinal) and L, M {L', M') refer to the quantum numbers of the total 
angular momentum of the initial (final) photon. Note that in Eq. (|27|) the dependence on 
the arguments q and q' also results from the dependence of Mj^ of Eq. (^Bp on p, given 
by P = —{q + if' )/2 in the center-of-mass frame. This implies that a naive substitution 
[p, ti;,g] [v, —uj\ —q'] in Eq. (P7| ) is no longer equivalent to photon crossing, which 
assumes that p is not affected by such a transformation. This can also be seen from 
energy conservation, where "naive" photon crossing would imply 

J = uo + E{-q) - E{-q') ^ -J + E{q') - E{q), 

which clearly contradicts the correct relation under crossing, u' —u. We will come 
back to this point in subsection III.C. 

Using the orthogonality property of the basis vectors V^(pLM, q), the multipoles can 
be extracted from Eq. (pT]) as 

j^l dq dq' V;{p'L'M\c^)M'^%Ms,u',q'-M,,u,q)V,{pLM,q). (28) 

Since the residual term is gauge invariant, q'^M^^ = = q^M'^ , it is sufficient to consider 
scalar, magnetic, and electric multipoles only |[30|1 , 

Mr(M^,^',g-";M„^,g) = 

E E 9p',W^{pL'M\q')H^t'''''''^\Mf,u:\\q'\-M,,u:M^^^^ (29) 

p,p'=0 LM 
L' ,M' 

where 

W^'ipLM, q) = VipLM, q) + 5po^V^{?,LM, q), p = 0, 1, 2. 

Finally, the dependence on the target spin projections is extracted by defining reduced 
multipoles,[] 

ffr'^")^(a;',|g'|;^,|g|)= E i-^)'^''^'^"'' ' Mj, I M,\S s) 

Mi,Mf M,M' 

x{L' M',L - M\S s)i/if™'''^*')(M^,^', \qlM,,u, \q\), 



(30) 



'^The right-hand side of Eq. (65) of Ref. contains a typographic error and should include a 
summation over the projections s |56||. 
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where / is the spin of the target, and any of the 25 + 1 projections s can be chosen. The 
selection rules due to the conservation of total angular momentum and parity are 



\L- L'\< S < L + L', 0<S<2I, 



-l)Vi+^'. 



From now on we assume the final-state photon to be real, u' = \q'\. The definition 
of the GPs requires to identify the leading-order behavior in \q'\ for any given multipole. 



In the Siegert limit of Ig*'! 
related by 



0, the final-state electric and longitudinal multipoles are 



lim H 



{2L',pL)S 
R 



lim 



'L' + 1 (3L',pL)S 



which in turn can be expressed in terms of -ff]^^ '''^^^ by gauge invariance. The leading- 
order behavior of a scalar or magnetic multipole of order L' is given by |^'|^ . The 
treatment of the initial-state photon is similar, except that one is interested in the behavior 
of the multipole for arbitrary values of The general relation between electric and 
longitudinal multipoles is given by 

8lM = ^Uq) + 



(31) 



(32) 



where £^j{q), £^^(g), and yij{q) are electric and longitudinal vectors of the multipole 
expansion, and vector spherical harmonics, respectively. While the first term on the right- 
hand side of Eq. (|32D can be expressed by the scalar multipole, the second term leads to 



the so-called mixed multipoles, which are neither of electric nor longitudinal type |30 



^r™^(M;,..',|g'|;M„.;,|g|) 
1 



(33) 



dg dg' V;{p'L'M\<^)M^^{Mf,u\q'-M,,u, \q\) [y^/+^(g) 

with reduced multipoles analogous to Eq. (|30D. The GPs of Ref. 0] are now defined as 

1 



p(p'L',pL)5^|^-.|) 



tj'^ \q 



Hi^'^'^'^'ico', \q\ 



uj'=0 



uj'^'\q 



(p'L',L)Sf , 



(P,P' = 0,1) 
(P =0,1). 



(34) 



Up to this point we have only incorporated constraints due to rotational symmetry, gauge 
invariance, and parity conservation. As such, the above analysis is valid in both a non- 
relativistic and a relativistic framework. However, additional restrictions apply, once a 
covariant, field-theoretical approach is chosen. In particular, if one assumes symmetry 
under particle crossing in combination with charge-conjugation invariance, one obtains 
a set of four independent linear equations, involving the ten generalized polarizabilities 

with v = 1 mm. 
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B. Expansion of Liu, Thomas, and Guichon 



Let us first reconsider the calculation of the GPs according to Ref. where the 
hadronic tensor M'^ is split into a leading and a recoil contribution. For that purpose 
one defines the matrix elements of the charge density and current operators of Eqs. (P^)- 
(0) between the ground state and the intermediate states X 7^ 0, 

Sx(g) = E^(^l^"e''''Ho). (35) 



a 



With this convention the current matrix elements of the direct-channel (subscript d) of 
Eq. (^ readg 



J,Mq) = {X\J{q,-qm = Px{q)+ttx{q) X q- ^ 

J%x{q) ^ {m-q')\x) = pW 



In the crossed channel (subscript c) the intermediate states propagate with momentum 



J^Mq') = (0|/(-r, -q')\X) = P*x{q') - ^^x{q') X q' - J^PW)- (36) 



Px = —q — q\ resulting in matrix elements depending on both q and q' . As in Ref. |^ 
we write these matrix elements as 

{x\j{-q\ -2q- q')\o) = iMq') + ^iMq, q'), 

{Q\J{q, -2q' - q)\X) = l^^xiq) + 5/c,ox(g, q'), (37) 

where 



4xoiq') 


= Pxi-q'] 




Jc,xoiq') 


= Pxi-r 


-iJ^xi-q") X q" 


Jloxiq) 


= Pxi-q) 




Jc,ox{q) 


= Pxi-q) 


— *■* 

+ zE^(-g) X q- 


"'^c,XO 


= SJlox = 


0, 


SJc,xo{q,q') 




-q'). 


SJc,ox{q,q") 




-q)- 



2M 



q 

2m' 



(3^ 



In the following we suppress the spin indices in our notation. 
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In Ref. [0, the leading part of the hadronic tensor is obtained by neglecting the 5Jc 
terms of Eqs. (|3^) in the crossed-channel matrix elements and by only taking account of 
the leading terms in the expansion of the denominators in a power series in oj' f\ 



-1 



Ef{-q>)+uj' -Exi'^) 
1 



X 



1 + 



00' 



AE 



X 



/2\ 



E 



x\ 



-q - q' 



UJ' 



Ei{q)-Ex{q) 



1 + 



1 



E,{q)-Ex{q, 



Ex{q. 



(39) 



where AEx is the excitation energy of the state X. As a result one finds [see Eq. (20) of 
Ref. [0] 



^R-LEADINci'l' ^ 9 , 



E 

Xy^O I 



Jiox{q')Jd,xM. 

AEx 



E 

Xj^O 



Ex{q)-E,{q) 



+ (40) 



The remaining contribution from the crossed channel is included in the recoil part, which 
also contains additional terms of first order in uj' from the expansion of the energy de- 
nominators in Eqs. (|55|). Accordingly, we obtain for the recoil part 



M^-RECOILiQ^Q. 



UJ' 



E / AZ7 \2'^dflx{^ )Jd,Xo{^) 

V 1 



Ex{q\ 



^cfiX 



-y ^ 

Exiq) - E,{q) 



Jcfix{q)^Jc,xM 



(41) 



This result differs from Eq. (21) of Ref. by the presence of the first sum on the right- 
hand side, the 5 J term of the second sum, and the 6J6J of the third sum of Eq. (^Tf) . 
However, this difference has no bearing on the calculation of the ten GPs with L' = 1. In 
this case we only need to analyze terms of Eq. (^) that are linear in c^;' = \q'\ [see also Eqs. 
(0)]. Furthermore, it follows from Eqs. ( pS] ) and (^) that the projection involves angular 
integrals of the type / dq'Y{'^.j,{q')M^ and / dq'Y{^,{q')M^^ for scalar and magnetic final 
dipole radiation, respectively. The recoil contribution of Eq. (^TJ) only modifies the GPs 
with a magnetic photon in the final state, as observed in Ref. [^. In fact, for = we 
need to consider J^qx, Jc,xoi ^^^d SJ'^xo- According to Eq. (^), the last term vanishes 



^In a nonrelativistic framework consistent with Eq. ([l|) one would have Ei{q 
—AEx- However, in Ref. [37| the relativistic expressions Ei{q) 
V'(M + l^ExY + q'^ have been used. 



Ex{q) 

v/Af2 + ^2 and Ex{q) 
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identically. The first two terms have to be always evaluated at q' = 0, because they are 
multiplied by expressions that are either manifestly or implicitly [SJcflx{Q,(l')] of order 
u'. From Eqs. (|36|) and (^) we see that the corresponding overlap integrals vanish due 
to the orthogonality of the excited states with respect to the ground state. In conclusion, 
the L' = 1 scalar GPs do not receive a contribution from Eq. (0). For the L' = 1 
magnetic GPs we need to consider fi ^ 0. In this case the terms proportional to u' 
do not contribute, because the angular integral of Y^^, multiplied by a g' -independent 
function vanishes. Furthermore, since SJc.oxiQ, q') is of order uj', the last term in Eq. (^) 
containing SJc,xo{(l,0) vanishes due to orthogonality. In conclusion, one is left with 

+ + (42) 
where the ellipses denote terms which do not contribute to the V = 1 GPs. 



C. Nonrelativistic 1/M expansion 

In this subsection we will consider a different procedure of ordering the contributions 
of the excited states to the Compton tensor. Our aim is to analyze the implications of 
photon-crossing symmetry for the GPs. We will see that some of the GPs will have to 
satisfy certain conditions for |g | — * 0. At first, we will not specify the reference frame but 
both photons will be allowed to be virtual [see Eq. (pT]) for X 7^ 0]. Only at the end we 
will assume the initial and final photons to be virtual and real, respectively, and restrict 
ourselves to the center-of-mass frame. 

The starting point is a separation of four-current matrix elements into an intrinsic 
current with respect to the center-of-mass system and a center-of-mass convection current 
53[| [see Eqs. (|1^-(|17D]. As an example, let us consider the following matrix element 



entering the direct channel 

(0|J^(-g^2p + g)|X) = (0|J^(-g',0)|X) + (0|5J^(-g',2p + g)|X), 

with analogous expressions for the other combinations. Here, we have defined for generic 
momenta a and h 

/i = : (0|5J°(a,6)|X) = 0, 

/i^O: (0|5J(a,6)|X) = A,(0|p(a)|X). 

Using nonrelativistic kinematics = p^/(2M), we separate the energy denominators 
into excitation pieces and kinetic contributions. 
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1 


1 


1 + 




-Ex{pf + q') 


^Ex 


AEx 


1 


1 


1 - 


UJ 


- Ex{pf - q) 




AEx 



I T{pf)-T{pf + q') ^ 



UJ 



/2 



AE 



X 



I TjPf) -Tjpf -q) ^ 



AE 



X 



i^Ex] 
(AEx) 



+ 0{3) 



+ 0{3) 



where C(3) refers to terms which are suppressed by l/{AEx)^, l/[{AEx)'^M], etc. In 
this representation we have 

Tfi.xi.^ A.P) = 

[(0|J^(-g',0)|X) + {0\5J^^{-q',2p + q)\X)][{X\r{q,0)\0) + {X\5riq,2p + q')\0)] 



E 

x^o 



X 



1 + 



UJ 



^ T{pf)-T{pf + q') ^ 



AEx 
{AEx] 



0{3) 



E 



AEx AEx 
mj%q,0)\X) + {0\5r{q,2p- q')\X)][{X\J^^i-q^ 



X 



UJ 



T{Pf) - T{pf - q\ 



AEx 

2 



UJ 



AE 



x 



AE 



X 



{AEx] 



0{3] 



(43) 



In order to expand Eq. ( P5D in powers of 1/M, we introduce the following abbreviations 
for the direct-channel {d) terms: 

A:^ = (X|r(g,0)|0), 
5A'', = {X\6r{q,2p+q')\0), 

B^, = {0\J^{-q',0]\X), 
SB^, = {0\6r{-q',2p + q)\X), 
AT, = T{pf)-T{pf + q'). 

The corresponding expressions for the the crossed channel (c) are obtained by simply 
replacing q —q' and fi ^ u. We note that both 6 and A terms count as order 1/M. 
The terms of leading order in 1/M and arbitrary order in 1/AEx, 



Tfr,x,Lo.i(l',<l) = - E 

Xj^O 



AE 



X 



uj' AEx + 



(44) 



do not depend on p. Equation (^4]) is symmetric under "true photon crossing," but also 
under "naive photon crossing" in any frame, because it is independent of p. As is shown 
in App. B, this property implies the following conditions on the leading behavior in 1/M 
for the reduced multipoles 



H 



(p'L',pL)S(_ y 
RI.o. 



{uj\ \q'\-uj, \q\) = Vp'^7pp(-)^+^'-^^il^^'^''(-^, l^'l). 



(45) 



Next we consider the 1/M corrections to Eq. (^4) to arbitrary order in 1/AEx- 
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E 



AEx - uj' 



+ 



AEx - uj' AEx 



AEx 

+ 



ATr 



UJ' 



AEx + ^ AEx + 00 



(46) 



The first line of Eq. (^) contains tlie modifications of tlie current matrix elements to first 
order in 1/M whereas the second line involves the center-of-mass kinetic energies of the 
final and the intermediate states, respectively. In this scheme, each order is separately 
photon-crossing symmetric. However, terms beyond leading order in 1/M will, in general, 
not be symmetric under naive photon crossing in the center-of-mass frame. 

We now turn to the case of real photons in the final state, uj' = \q'\, and explicitly 
specify the center-of-mass frame. As we are not interested in terms beyond linear order 
in cj' and 1/M, we make use of energy conservation. 



UJ u 



q- 



2M' 



to expand the crossed-channel terms as 
1 1 



AEx + 00 



X 



LO 



+ 



io'q^ 



AEx 2MAEx {AExfM 



l2\ 



We then obtain from Eq. (0) up to and including terms linear in uo' and 1/M 



Tfi;X,lA^'^^) 



E 



AEx 

2 



1 + 



UO 



AE 



X , 



B^,A^, + 



uo 



R'" 4'' 

2M{AExY " ' 



uo'q'^ 



(AEx^M 



AEx 

b':a^ 



AE 



X , 



b:a^. 



(47) 



Note that via energy conservation, Eq. ( ^41) has also generated terms in Eq. (^Tf ) which 
are of higher order in 1/M. 

For further evaluation of Eq. (|46|) we use 



and obtain 



'^fi;X,n.l.o.i? ^Q, 



AT, = 0{u;'^), AT, 



1 



E 

Xy^O 
1 

AE^ \ 

._i Ri' /IM 



AEx 



2M 



UJ 



AE 



X 



UJ 



M 



[B^,5A^, + 5B^,A2 



AEx, 
1 

{AExT^ 



[B^jA^^ + 6b::a^j 



UJ 



AE 



X 



q-q' 



M ' "{AExY 



B^At- 



+O{uo", 1/M 



(48) 
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The last two terms of Eq. (|47D cancel with two corresponding terms in Eq. ( ^Sf) and the 
sum can be written as 



E 



AE 



+ 



X 



AE 



X 



2'^d^d 



+ E 



1 



M 



-E 

XytO 

-E 

x^o 



AE 



+ 



B^At 



B':6At^ + 6B':Aii 



X 



AE 



X 



OJ' 



{Bi;SA^, + 5B^,A:, - B:5A^^ - 6B'^A>i) 



+O{uo'\l/M'). 



(49) 



As in the previous case we will now analyze which terms of Eq. ( ^Q]) actually generate 
a contribution to the L' = 1 GPs after angular integration with the spherical harmonics 
^im'(^')- Again, all terms explicitly proportional to uj' will not contribute, because at 
leading order they are multiplied by expressions which do not depend on the direction of 
q' . The term vanishes for = and is of order a;'^ for /i 7^ 0, and one is left with 



'^fi;x{?'^ Q) — '^fi;X,LEADING{^' ^ ^) + ^/i;^ 



fi;X,RECOIL 



where 



fi;X,LEADING\1 ; 1 . 



'^fi;X,RECOILi'l ; ^, 



-E 

x^o 

-E 
+ E 



B^ B^ 
AEx AEx 
B^JA^, + 5B:A^^ 
AEx 
1 q-q' 



(50) 



(51) 



x^oi^Exr M 



B:A>i 



E 



B'iSAit 



x^o AEx 

+■■■+o(c^'^l/M2), 



(52) 



where the ellipses refer to terms which do not contribute to the V = 1 GPs. The terms in 
Eq. (|5lD generate contributions to the GPs which result entirely from intrinsic currents. 
The first line of the recoil term in Eq. (|52|) corresponds to 1/M corrections of the virtual- 
photon absorption vertices in the direct and crossed channel, respectively. The second and 
third lines are 1/M corrections of the crossed-channel energy denominator and real-photon 
vertex, respectively. These last two corrections only affect GPs involving a magnetic 
photon in the final state. 
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Finally, by taking the nonrelativistic limit of the energies in the scheme of Ref. |37 
Ei{q) = Jm^ + g 2 ~ M + ^ 



2M' 

2 



Ex{q) = = J Ml + ~ M + AE^ + ^ 



2M' 



it is straightforward to verify that the two expansion schemes coincide up to and including 
terms of order uj' and 1/M. They only differ in their separation into leading and recoil 
terms. 



IV. GENERALIZED POLARIZABILITIES IN A NRCQM 



In this section we discuss the GPs of the nucleon in the framework of a nonrelativistic 



system consisting of three constituent quarks. As in Ref. we restrict ourselves to the 
inclusion of the A(1232) resonance and the low-lying negative-parity baryons Di3(1520), 
Sii(1535), S3i(1620), Sii(1650), Si3(1700), and D33(1700). 



A. Matrix elements in the Isgur-Karl model 



To be specific, we employ the model of Isgur and Karl ||5^ which describes the quark- 



quark potential by a harmonic-oscillator term plus a spin-dependent hyperfine interaction. 
The baryon states are expressed in a basis of SU(6) harmonic-oscillator wave functions, 
with the SU(6) multiplets generated by the combination of SU(2)spj„ and SU(3)//a„or 
multiplets. In particular, the nucleon and the A(1232) resonance belong to the ground- 
state spin-1/2 octet, ^8, and the spin-3/2 decuplet, ''lO, of the 56 SU(6) supermultiplet, 
respectively. The multiplet of states associated with the negative-parity orbital excitation 
is classified in terms of a 70 supermultiplet of SU(6) which decomposes into ^1, ^8, ^8, and 
2 10 multiplets. According to Ref. |5^, the strength of the hyperfine interaction is fixed to 



reproduce the experimental mass splitting of and A(1232) states, while the remaining 
orbital excitations of the 70 multiplet are constructed with mixing parameters describing 
the empirical spectrum quite well. 

Since we are interested in the non-strange sector only, we assume that all three quarks 
have equal masses m^. In addition, when calculating the matrix elements of the electro- 
magnetic current of Eqs. ([161) and (p^Tj), we take advantage of the overall symmetry of the 
SU(6) harmonic-oscillator wave function. This allows us to simplify the matrix elements 
of one-body operators, 

{A\j2d,\B) = ^A\d,\B). (53) 
1=1 

As a result, the overlap integrals of Eq. (|35D can be written as 

px{q) = 3 / dpdXe-'^^'^-^ (P^x Q <PN, (54) 
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ZrUq J 

where 

4>N = <^7v(p, A, ^,M, ^,r^), 

0x = (pxip, A, Jx, Mx, /x, Tx) 

denote the internal wave functions of the nucleon and the excited states, respectively, 
with an obvious notation for spin and isospin labels. We have introduced the standard 
Jacobi coordinates p = (n — r2)/V^ and A = (ri + r2 — 2r3)/^/6, and made use of 
r 3 = rs — ^ = — -y/2/3 A. Furthermore, Q = {1/6 + and (T3 denote the charge 

operator and the Pauli matrices of the third quark, respectively. 

Explicit expressions for the contributions of the P-wave negative-parity states are 
given in Appendix 0. 



(55) 
(56) 



B. GPs in the framework of Liu, Thomas, and Guichon 



We first discuss the results for the proton GPs that we obtain with the same con- 
ventions used by Ref. [^| for the separation into leading and recoil terms. Following 
Ref. we use rriq = 350 MeV for the quark mass and a = 320 MeV for the oscillator 
parameter. As was pointed out in Ref. |37], the proton polarizabilities do not receive any 
contribution from the A^(^8) multiplets. The mixing parameters ax encoding the A^(^8) 
and A(^IO) composition of the resonant states are taken from Ref. and are listed in 
Tab. I. 

With these assumptions we find for the leading contributions to the Compton tensor 
of Eq. 



S,Jx 




^2,S,Jx 



E{q)-Ex{q), 



para 



,(11,11)5 



(57) 
(58) 

(59) 
(60) 



The angular coefficients Z of the leading contributions are given in Tab. II. The diamag- 
netic term of Eq. (^) originates from the modified seagull term S^^ and contributes in 
the spin-independent case only. The mixed GPs are given as the sum of two terms, cor- 
responding to the contributions from the convective (Pp'^'^^'^) and spin [P^^'^^^) terms of 
the current at the virtual-photon vertex 
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p(01,l)5 ^ p(01,l)S p{01,l)S 



,(01,1)5 



,(01,1)5 



108^3 m,a2 



36V3 rriga^ 



E 



X 



+ 



M-Mx E{q)-Ex{q), 



X=N*,A* 



2L- 



7l,5,Jx 



+ 



2'l,5,Jx 



M-Mx E{q)-Ex{q) 



(61) 
(62) 

(63) 



In the scheme of Ref. ^ the spin-dependent GPs p(ii'02)i^ p(ii,oo)i^ p(ii,2)i^ all of 
which lead to Ml radiation in the final state, vanish identically at leading order. 

The recoil corrections are exclusively generated by the crossed-channel diagrams and 
only modify the GPs with a magnetic final photon (see the discussion at the end of 
Sec. III.B). To start with, the two GPs p("'00)i and p(ii'02)i 

receive a non- vanishing recoil 

contribution 

1 



P. 



(11,00)1 



recoil 



(11,02)1 



recoil 



a/3 



1 1 
^3 



E 



"X^llOO 



Ex{q)[E{q) - Ex{q)f 



1 - 



Ex{q)\E{q)-Ex{q)\ 



3a2 



(64) 



„2 yJx 
"'X^1102 



X=N 



i^^Ex{q)[E{q)-Ex{q)f [ 



Ex{q)[E{q) - Exjq)] 
3a2 



(65) 

When discussing the recoil contribution to the remaining three polarizabilities P^HolP^ 
and Pllloii^ 1 it is useful to distinguish between terms which result from the spin- 
independent (C) and the spin- dependent (S) part of the virtual-photon vertex. We find 



3(11,11)5 
recoil 



(11,11)5 



recoiLC 



3(11,11)5 
recoiLT, 



p(ll,ll)5 p(ll,ll)S 
recoil,C ' recoil,T, 5 

18 m2 



(66) 



2 ^S,Jx 
"X^llll 



18 rrig 



-q^ /Ga 



X=N*A* 

Ex{q)[E{q) - Ex{q)f 



2Ex{q)[E{q) - Ex{q)] 
3a2 



(67) 



E 



a 



S,.Jx 
1111, s 



X' 



x=N'A* 24 Ex{q)[E{q)-Ex{qW 

Ex{q)[E{q) - Ex{q)] 
3^2 



(68) 



P, 



i(ll,2)S 
recoil 

(11,2)1 
recoil, C 



P, 



(11,2)1 



recoil, T, 



p(ll,2)5 p(ll,2)5 

recoil,C ' recoil,T,^ 

1 1 



6^5 m2 



E 



„2 yJx 
"X"^1102 



(69) 



£;x(g)[£;(g)-^x(g)]^ 



Ex{q)[E{q) - Ex{q)\ 



3a2 



1 1 



L102 



2v^rn2 ^ 2/. Ex(g)[i^;(g)-i?x(g)]^ 



Ex{q)[E{q) - Ex{q)] 
3a2 



(70) 



(71) 
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The values for the angular coefficients Z of the recoil contributions are given in Tab. III. 



C. Comparison with Liu, Thomas, and Guichon 

We now compare our results with those of Ref. As has been stated in Sec. III. A, 

the spin-dependent polarizabilities differ by an overall factor of three due to the different 
definition of the reduced multipoles in Eq. (|30|). In contrast with Ref. [0, we find, at 



leading order, two simple relations for the angular coefficients of the direct and crossed 
channels [see Eqs. (p7D-(|6^)], 



^f/'"^ = (-l)^^'^a/''"- (73) 

In particular, as a result of Eq. ( fT^j ) combined with the relative phases between the direct- 
and crossed-channel contributions to P^f^a^^"^ ■• leading term of the p(o^'°^)^ and p(^^'^^)^ 
polarizabilities vanishes at the real-photon point. This is in agreement with the constraint 
by photon crossing as derived in Eq. (p3|) of App. B, but in contrast to the results of 



Ref. In addition, in our calculation the diamagnetic contribution -Pjj^a'^^^^ is only 



3/7 of the result of Ref. P^. We also find different expressions for p(oi>i)'S'^ particular 
for the relative phase of the direct- and crossed-channel contributions. Furthermore, the 
angular coefficients 2'/|oo -^no2 occurring in P^^f^^a'^ ^'^d pHIm^^ smaller by a 
factor of 1/2. In addition, we find a recoil contribution to PrlcolP^ both spin-flip and 
no-spin-fiip transitions, while in Ref. such a contribution is absent for S = 1. 

The numerical results of Eqs. (^7|)-([7l|) are shown in Figs. 1 and 2 together with the 
calculation of Ref. p^. For the generalized electric polarizability 



the two results are in agreement (full line of Fig. 1). The discrepancies in p(oi'Oi)i^ p{ii,ii)i^ 
and p(oi'i2)i 

originate from the contributions of the crossed channel relative to the direct 
one. The different results for the generalized magnetic polarizability 



3 



are mainly due to the discrepancy in the calculation of the diamagnetic term, while 
the different evaluations of the recoil terms result only in small deviations. Finally, the 
differences in calculating the p(oi'i)'5' and Prll^]^ polarizabilities give rise to discrepancies 
of almost one order of magnitude. In particular, the P^ll^]^ polarizability receives its 
main contribution from -Precoi/s which has been neglected in Ref. 



^'^We stress that both calculations start from the same model and use the same parameters and 
approximations . 
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D. Comparison with other calculations and experiment 



In Figs. 3 and 4 we compare the two expansion schemes of Sees. III.B and III.C. 
In each graph the sohd line represents the full result according to the scheme of Liu, 
Thomas, and Guichon. Recall that in this framework the energy denominators of the 
crossed channel contributions are written using relativistic kinematics. The dashed line 
corresponds to a consistent nonrelativistic expansion up to and including terms linear in 
1/M. The relevant expressions can be found in App. 0. The contributions of leading 
order in 1/M [see Eqs. (|5T| ) and (^2])] are separately displayed as dotted lines. 

First of all, we note that due to photon-crossing symmetry the leading contributions to 
both p(oi'Oi)i and p(ii'ii)i vanish identically. The use of relativistic expressions in the en- 
ergy denominators for the crossed-channel terms leads to pronounced differences between 
the two schemes, as soon as the leading term is vanishing or small. A striking example 
is given by the difference between the solid and dashed lines in p(oi'Oi)i of Fig. 3, which 
is entirely due to this different treatment of the crossed-channel energy denominators. 
Finally, at the real-photon point, the leading contributions of both expansion schemes are 
equal, whereas the recoil terms differ by the contribution of second-order terms in 1/M. 

Figures 5 and 6 display our results (full lines) together with the predictions of the 
linear sigma model (dashed lines), an effective Lagrangian model ||39| (dotted lines), 
and heavy-baryon chiral perturbation theory (dashed-dotted lines). 

An unpolarized measurement can be analyzed in terms of three structure functions 
Pll, Ptt, and Plt HS'H which are products of the GPs and the electromagnetic 
Sachs form factors Ge and 6**^,0 



PLLiWl) 
PTTiWl) 

PMlql) 



-2V6MG^(Qg)p(°i'°i)°(|g|), 



(74) 



'-G,M) {2u;oP^'''''^\\q\) + V2\qnV^P^''''^\\q\) + P^^^'^^^^dg |)]} 

3GM{Ql)\mV2P^'''''^\\q\)-P^'''''^\\q\)/coo], (75) 
[3M\q 



GE{Ql)P^'''''^Mq\) 



+—^Gm{qI) 




p(n,oo)i^|^-.|) + ^p(n,02)i^|^-.|) 



GE{Ql)P^'''''^'{\q\) + l^^^^GM{Ql)P^'''''^\\q 



(76) 



where uq = i-i-'|a;'=o = M — yM^ + |g ^ and Ql = —q'^\^>^Q = —2Mujq. We note that 
the second equations of Eqs. ([75| ) and (|76D , respectively, rely on symmetry under particle 



crossing and charge conjugation [see Eqs. (21c) and (21d) of Ref. ||3^] which is not satisfied 



^^The fourth structure function P'^j. of Ref. |3C| is related to P^t if symmetry under particle 
crossing and charge conjugation is applied |]35|. 
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in the NRCQM. Table IV contains the predictions for the three response functions for |g | = 
(0,240,600) MeV, corresponding to the real-photon point, MIT-Bates and MAMl 
|^6[| kinematics, respectively. For the Sachs form factors we used the parametrization 
of Ref. |]58|. For Ptt and Plt we quote both results obtained from Eqs. ( |75| ) and (|76D, 
respectively. 

Finally, in table V we compare the predictions with the first experimental information 
obtained at MAMl for the linear combination Pn — Ptt/^ and Plt at Ig*! = 600 
MeV. 



E. Particle crossing, charge conjugation, and gauge invariance 



The original definition of the GPs of Ref. |30] was based on angular momentum con- 



servation, parity conservation, and gauge invariance. In Refs. ||3^J35[] it was shown that 
only six of the originally ten GPs are independent if particle-crossing symmetry in com- 
bination with charge- conjugation invariance is imposed. As a result a set of four linear 
equations was obtained. 



p(01,l)0 



P 



(11,2)1 



/ 2p(oi,oi)o _^ J_p(ii,ii)o 
3 \/Q 



^p(ll,ll)l + ^3^^^(11,02)1 



p{01,l)l = _ 1 .^cUoP^^^'''^' + ./^^P(ll'll)l - |2p(01,12)l 

|g r V 3 V 3 cuo V3 



p(ll,00)l ^ ^73 1^1 p(01,01)l L|g'|2p(ll,02)l 



(77) 



(7J 



(79) 



^0) 



where Uq = uj\^/=o. Following common use, we choose for the six independent GPs a, /?, 

p(01,01)l^ p(ll,ll)l^ p(01,12)l^ p(ll,02)l gi^g^ pigg^ 3^ 5_ 

However, in a nonrelativistic framework particle crossing is not a symmetry of the 
Compton tensor (see Sec. 2.3 of Ref. |§5[)- Thus one cannot expect the relations of 
Eqs. (^ - (Hg) to be satisfied in the NRCQM. Having defined {i = 1,2,3,4) as the 
difference between the left-hand and right-hand sides of Eqs. ( [77| ) - ( |80| ) normalized to 
their right-hand sides, we show the discrepancies Aj in Fig. 7 as function of |g |. Clearly, 
the relations of Eqs. (^) - ( ISDD are strongly violated on the average. 

Finally, another important limitation of the NRCQM is due to a violation of gauge in- 
variance. This problem can be traced back to essentially two causes: first, the Isgur-Karl 
model includes some effects of the anharmonic terms in the qq potential only perturba- 
tively in the energy but not in the wave functions. Such a treatment leads to a mismatch 
between the resonance masses in the energy denominators of the Compton amplitude 
and the baryon states which enter in the current matrix elements. Second, the actual 
calculation truncates the configuration space to only few intermediate states, while gauge 
invariance requires, in principle, the full set of intermediate states. The well-known exam- 
ple in this context is, of course, the low-energy Thomson limit which requires the inclusion 
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of the sum over all electric-dipole excitations (including the negative-energy states of a 
relativistic theory). 



V. SUMMARY AND CONCLUSIONS 



We discussed the general form of the virtual Compton scattering tensor for a nonrel- 
ativistic composite system. In particular, we focussed the attention on the generalized 
polarizabilities of the proton, defined from the multipole expansion of the Compton tensor. 
We performed a consistent nonrelativistic expansion of the structure-dependent amplitude 
allowing us to identify the constraints due to photon crossing. As a model calculation, 
we reconsidered the proton GPs in a nonrelativistic constituent quark model. The model 
satisfies the constraint due to photon crossing at the real photon point, but does not 
provide the relations among the GPs due to nucleon crossing in combination with charge 
conjugation. As a consequence of its limitations regarding relativity, gauge invariance and 
chiral symmetry, the results of the model should be treated with some care. There clearly 
is room for improvement in any of the above-mentioned shortcomings. Nonetheless, the 
predictions provide an order-of-magnitude estimate for the nucleon resonance contribu- 
tions and as such are complementary to the results of the linear sigma model and chiral 
perturbation theory emphasizing pionic degrees of freedom and chiral symmetry. 
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APPENDIX A: GAUGE-INVARIANT MODIFIED GROUNDSTATE POLE 

TERMS 

Here we derive the result of Eq. (^) for G^*^ which, in combination with the ground- 
state pole terms, constitutes the gauge- invariant tensor Mp^ . We start with the expression 



[see Eqs. (8a, b) of Ref. with the replacement k ^ q' and k' —q] 



Tf-{q\q^p) = -t{2nfj dhe''''-^{QMfPf\T[J^iz)r{Q)]\OMS): 

where J^{x) = exp{iHQt)J^{x) exp{—iHQt) . In general, we allow for a change in the spin 
projection from Mj to Mf. The groundstate pole contribution is given by 

-i{27TfJ2 I I d^ze''^'" \0{zo){OMfPf\Jf'{z)\OMP){OMP\r{0)\OMS) 



M 



+9{-zo){OMfPf\r{0)\OMp){OMp\J^'{z)\OMS) 



(Al) 
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The following procedure is very similar to the one used in deriving Ward-Fradkin- 
Takahashi identities [p^-pT| in Quantum Field Theory [see, e.g., Chap. 6.1 of Ref. p2| j. 
Let us contract Eq. ( [All ) with q'^ (arguments suppressed), 

Q',Tj^:f^ = -^(2^)' Hld'P I d'z [-^^,e^^'-) [■■■]. 

Symbolically this expression is of the type 



where we made use of a partial integration, and assumed that the interaction is "adiabat- 
ically" switched on and off to get rid of the surface terms at zq = ±00. Similarly, use of 
the divergence theorem has been made. After applying this result to the above case. 



M 

we use the relation 



Q',nj,= i2nrT.Jd'Pjd'ze^'^'-d,[. 



d,[9{±zo)g''{z)] = ±5izo)g'{z)+9{±zo)d,g'^iz), 

and obtain 

g;ri^}, = (27r)=^^ / d^P j d^ze"^'-' [e{zo){OMfPf\d^J^'{z)\OMP){OMP\r{0)\OMS) 

+e{-z^){mfPf\rmmp){mp\d^j^{z)\m,p,) 

+6{zq) {OMfPf\p{z)\OMP) {OMP\ r{0)\OMS) 
-6{zo) {OMfPf I r{0)\OMP) {OMP\p{z)\OMS) • 

With d^J^{z) = as an operator identity, the first two terms on the right-hand side of 
the equation vanish. Performing the integration with respect to 2:05 ^'Pplying translational 

invariance as p{z) = exp{—iP ■ z)p{0) exp{iP ■ z), integrating first with respect to z and 
then with respect to P, the two remaining terms yield 



q'^T^f^ = {2nfY.[{^MfPf\p{0)\0Mpf + q'){OMpf + q'\riO)\OMS) 

M 

-{OMfPf\r{0)\OMp,, - q'){OMp,, - q'\p{0)\OMS)] • 

Finally, the operator p of Eq. (^) is diagonal in the spin projections and we obtain 

g^Ti^}, = (27r)^ [(OM;p^|p(0)|OM;p; + q')mjpf + q'\rmOM,p,) 

-{QMfPf\r{0)\OMS - q'){OMS - q'\p{Q)\QMS)] • (A2) 

First we consider Eq. ( [A2| ) for u = 0: 

Q'.T^fi = [(OM;|p(-g')|OM;)(OM;|p(g)|OM,) - (OM;|p(g)|OM,)(OM,|p(-g')|OM,)] = 0, 

(A3) 
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where we inserted Eqs. (|T4|) and ([T7|) and made use of the diagonal nature of p. Similarly, 
inserting Eqs. (0) - ( p!7D for u = j we obtain 



(OM,|p(-g')|OM,) 
= {OMf\r'^{q)\OMi) [{OMM-Q')mf) - (OM,|p(-g-")|OM,)] 
+^{0Mf\pi-Q')mf){0MMQ)mf)5M,M, 

= ^{0Mf\pi-Q')mf){0MMq)m,)6M.M,: (A4) 

where we made use of the fact that the groundstate matrix elements of p are diagonal 
and do not depend on the projection. 

The calculation of q^Tpji proceeds in a completely analogous fashion. Equations (|A3| ) 
and (|A4|) suggest to add the term of Eq. (^51) , 



G^'{q',q) = G'^{q\q) = 0, G'^{q',q) = 5.,^(0M;|p(-g')|0M;)(0M,|p(g)|0M,)5M,M,, 

with the result that Mp^ = Tp^j^ + G^'^ is gauge invariant. In particular, Mp' depends on 
groundstate properties only. 



APPENDIX B: PHOTON CROSSING CONSTRAINTS AT LEADING ORDER 

IN 1/M 

In Sec. III.C we have seen that the leading term of the residual amplitude, in a 1/M 
expansion, satisfies the photon crossing constraint [see Eq. (^41)] 

2^/^k.o.(?',?)=2^;k.o.(-9,-?')- (Bi) 

Because Eq. does not depend on the average target momentum, it is also symmetric 
under "naive photon crossing. 
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T£,,.(M;, u', q'- M,, q) = T^^Mf, -co, -q; M„ -u', -q'), 
which in terms of the multipole expansion of Eq. ( p7D implies 

4vr E gp'p'V^ip'L'M',r)H^/^iT''''''''\Mf,u;', \q'\;M.,u, \q\)V'*ipLM,q)g^ 



p,L,M, 
p',L',M' 



PP 



4vr E appVipLM, -q)HS.y \Mj, -u;, |g |; M„ -uj\ \q'\)V^*{p' L' M\ -q')9,',' 

p,L,M, 
p' ,L' ,M' 



4vr E {-f^''^''^'''V^{p'L' - M',q')H'^iy'''''\Mf, -u, |g |; M„ -u' , \q'\) 

M, 
,M' 

xV''*{pL-M,q)gp,p. 



p,L,M, 
p',L',M' 
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In the last step we made use of 

V^^*{pLM, -q) = ^7pp(-)^+^K^(pL - M, g). 

With the orthogonahty property of the muhipole basis {V^{pLM, q)} (see appendix C of 
Ref. 13) we find 

(-)^^''+"'-'''^7pV'^?Pp4t.*''''^'"'''HM/, |g I; M„ (B2) 

Because of the orthogonahty relations of Clebsch-Gordan coefficients we finally obtain for 
the reduced multipoles 

H^/Z:''^'>'{Mf,u',\q'\-,M.,uj,\q\) = 

(-)'^+^"Vp'^?pp^S'^'^''(M/, -u, \q\; M„ -u', \q'\). (B3) 

We stress that the above derivation holds only for those pieces that are independent of p. 
It is only in this case that true photon crossing is equivalent to naive photon crossing in 
the center-of-mass frame. 



APPENDIX C: OVERLAP INTEGRALS OF THE CURRENT OPERATOR IN 

THE NRCQM 



The contribution of the SU(3) multiplets to the overlap integrals defined in Eqs. (|5 
^ ) is given by (g = |g |) 

(Cl) 



(C2) 

2 



1 h Jx 



S = i^^le~^'/^"Vl + 4r l Vr-n^+^^^A/2Jy + l ( ^ 

xnw-p('?)(xii^3ixt,), (C3) 



3 1 

2 



xnW-p(9)(x!l^3|xiJ, (C4) 
Ejg-) = e-V6.^(xija3|xtj, (C5) 
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1 + 4r^ '^t^^) 



(C7) 
(C8) 



where Cm is the spherical basis vector and a the oscillator parameter. The eigenstates 
of the total spin S of the three quarks have been denoted by x^J"^ for S* = 3/2, and 
for S = 1/2, with A indicating symmetry under interchange of the (12) quark pair. We 
note that the definitions for the overlap integrals introduced in Eqs. (p4D-(^) differ by a 
factor 3 from the corresponding expressions in Eqs. (9)-(ll) of Ref. |3^. In addition, we 



found a different result in the explicit calculation of the integral entering into S^^2g)(^) 



(the integral in Eq. (14) of Ref. is smaller than our result by a factor 2-\/2), while we 
agree with the results for the integrals contributing to the remaining terms given in Eqs. 



(12) and (13) and Eqs. (15)-(19) of Ref. |37 . 



APPENDIX D: POLARIZABILITIES IN THE NONRELATIVISTIC l/M 

EXPANSION 

In this appendix we collect the results for the polarizabilities obtained from the mul- 
tipole expansion of the Compton tensor in the l/M nonrelativistic limit. 

The leading contributions corresponding to the terms in Eq. (|5TD read (g = |g |) 

V/6a^ y ^2 L^i fad ^2) 

a + (-i)^), (D3) 



p(01,01)5 ^ 


1 1 -g2/6a2 

18 a2 


p(01,12)l ^ 


1 /s 1 _ 
18 V 5m,a2^ 


p(ll,ll)5 _ 
para 




p(ll,ll)S _ 
dia 


. 1 1 
3V6 "^<7« 


p(01,l)S 


1 1 


36^3 m^a^ 



(D4) 



Comparing these results with Eqs. (|F7|)-(|52D for the leading contributions of the po- 
larizabilities in the convention of Ref. fS^], we notice that pCoi'Oi)^, P^lla^^^ , and p{oi'i2)i 



involve the same matrix elements of the current in both expansion schemes, while for the 
propagator in the crossed channel we now have (M - Mx)~^ instead of [E{q) - Ex\~'^. 

The diamagnetic contribution to p(ii'ii)o from the two-photon interaction does not 
change in the l/M expansion. 

The polarizability pC^^'^)-^ has a leading contribution only from the spin-dependent 
term of the current, p|°^'^'''^, while the contribution from the convective current corre- 
sponds to l/M corrections that are taken into account in the recoil term. 

The recoil contributions can be obtained from the respective terms of Eq. (^), 
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5(01,1)5 
recoil 



3(11,00)1 
recoil 



P. 



(11,02)1 



recoil 



p(ll,ll)5 
recoil, C 



3(11,11)5 
recoil, 'E 



p(ll,2)l 
recoil, C 



f>(11.2)l 
recoil, 'E 



V2 1 

108^3 mga2' 



E 



X=N*A* 



'yS,,Jx 

M-M 



X 



1 g2 



a/3 
1 1 

V3 "^g 



-cp- l%c? 



E 



"x^iioo 



M[M -M, 



X 



E 



"X^1102 



E 



"X^1111,F 



E «- 

M[M -Mx 



^, M[M-MxY 



1 - 



3a2 
M[M-Mx] 

2M[M-Mx] 



3a2 



aiii,s 
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M[M 



1 1 



6\/5 m2 



3a2 



«2 yJx 

"'X^1102 



Mxl2' 



2^5 



E 



2 



-11 



^1/2 



^1102 



24 M[M-M 



X 



3a2 



(D6) 
(D7) 
(D8) 
(D9) 



(DIO) 
(Dll) 



(D12) 



Comparing with the results in Eqs. (|6^-(|7TD, we find the same expressions 

fQj. p(ll,00)l 

p(ii,02)i^ and P("'ii)i after the substitutions [E{q) - ^x(g)]"^ [M - Mx]'^ and 
Ex\q)[E{q) - Ex{q)]-^ M-'[M - Mx]-'. 

As noticed previously, the additional recoil term for p(*^^'^)i corresponds to the p^^'^''^ 
term in Eq. (0). Furthermore, Prll^]]j does not contain the contribution from the 
convective current at the virtual-photon vertex resulting from the term proportional to 
'^c,fx{Q)^Jc,xi{^') Eq. (|4^), because this contribution corresponds to a higher-order 
correction in 1/M. 
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TABLES 

TABLE I. Mixing parameters for the ^8 and ^10 representation of the P-wave baryon reso- 
nances. 



X 


A(i ,1535) 


,1650) 


A(| ,1520) 


iV(| ,1700) 


A(i ,1620) 


A(| ,1700) 


ax 


0.85 


0.53 


0.99 


0.11 


1.0 


1.0 



TABLE IL Angular coefficients for the leading contribution to the GPs. 



L 


S 


Jx 




ZS,Jx 


yL,S,Jx 


vL,S,Jx 
'-'ac 


zi 


1 





1/2 


V273 




-2/V3 


2/V3 




1 





3/2 




2^/2/3 


2/V3 


-2/V3 




1 


1 


1/2 


2/3 


-2/3 


-2V2/3 


-2^/2/3 




1 


1 


3/2 


-2/3 


2/3 


-V2/3 


-V2/3 


-1 


2 


1 


3/2 






V30/3 


\/30/3 





TABLE in. Angular coefficients for the recoil contribution to the GPs. 



s 


Jx 


yJx 
^1100 


^1102 


r^StJx 
^1111, C 


^1111, E 





1/2 








IV6 





3/2 






Vi 


-1^6 


1 


1/2 


1/27 


l/27\/2 


-1/3 


-1/6 


1 


3/2 


-1/27 


-l/27\/2 


1/3 


1/6 



TABLE IV. Structure functions Pll, Ptt-, and Plt in GeV ^. The two entries for Ptt and 
Plt originate from the first and second equations of Eqs. ([7^) and ([76|), respectively. 





Pll 


Ptt 


Plt 


\q \ = MeV 


37.0 


-0.1 


0.0 


-11.2 


- 15.8 


\q \ = 240 MeV 


28.7 


-2.8 


- 1.4 


-8.8 


- 12.4 


\q \ = 600 MeV 


9.9 


-5.8 


- 3.1 


-3.2 


-4.5 
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TABLE V. Structure functions PLL-PTr/e and Plt for |^| = 600 MeV in GeV-^ (e = 0.62). 
The two entries for P^l — Ptt/^ and Plt in our calculation originate from the first and second 
equation of Eqs. (75) and (^), respectively. 





Pll - Prr/e 


Plt 


This calculation 


19.2 14.9 


-3.2 -4.5 


NRCQM of 13 


11.1 


-3.5 


Experiment 46] 


23.7 

±2.2 ±0.6 ±4.3 


-5.0 
±0.8 ±1.1 ±1.4 
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FIGURES 




FIG. 1. GPs in the NRCQM as function of virtual-photon momentum q = \q\. Full lines: 
our results in the scheme of LTG, taking account of leading and recoil terms; dashed lines: 
calculation of Ref. |37]. Note that the two calculations coincide in the case of a. 
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FIG. 2. GPs in the NRCQM as function of virtual-photon momentum q = \q\. Full lines: 
our results in the scheme of LTG, taking account of leading and recoil terms; dashed lines: 
calculation of Ref. |37]. 



33 




FIG. 3. GPs in the NRCQM as function of virtual-photon momentum q = \q\. Full lines: 
our results in the scheme of LTG, including the leading and recoil contributions; dashed lines: 
nonrelativistic 1/M expansion of the Compton tensor, obtained from the sum of the leading 
and recoil terms; dotted lines: leading contributions. Note that dotted lines are only included 
if the leading contribution is finite and different from the total result. To the order considered, 
the polarizabilities a, p(oi:i2)i^ ^^^^ p(oi,oi)i receive only contributions from the leading terms 
of Eq. (51), whereas p(ii'02)i consists only of a recoil contribution from Eq. (52). 
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FIG. 4. GPs in the NRCQM as function of virtual-photon momentum q = \q\. Full lines: 
our results in the scheme of LTG, including the leading and recoil contribution; dashed lines: 
nonrelativistic 1/M expansion of the Compton tensor, obtained from the sum of leading and 
recoil terms; dotted lines: leading contributions. Note that dotted lines are only included if the 
leading contribution is finite and different from the total result. To the order considered, 
receives only a contribution from the leading terms of Eq. (51), whereas p(ii'00)i g^j^j p{ii,2)i 
consist only of a recoil contribution from Eq. (52). 
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FIG. 5. Results for GPs in different model calculations as function of squared momentum 
transfer at = 0, Qq = Q'^\uj'=o- Full lines: our results in the NRCQM with the scheme 



of LTG; dashed lines: linear sigma model [Q; dotted lines: effective Lagrangian model |39|| ; 
dashed-dotted lines : heavy-baryon chiral perturbation theory E 
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FIG. 6. Results for GPs in different model calculations as function of squared momentum 
transfer at oj' = Qq = Q'^\uj'=o- Full lines: our results in the NRCQM with the scheme of LTG; 
dashed lines: linear sigma model [EH]; dotted lines: effective Lagrangian model [El]; dashed-dotted 



lines : heavy-baryon chiral perturbation theory |44]. 
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FIG. 7. Violation of the relations of Eqs. ( |77[) - (|80D as function of q 
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